Cutoff-independent regularization of four-fermion interactions for color 

superconductivity 
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We implement a cutoff-independent regularization of four-fermion interactions to calculate the 
color-superconducting gap parameter in quark matter. The traditional cutoff regularization has 
difficulties for chemical potentials n of the order of the cutoff A, predicting in particular a vanish- 
ing gap at /i ~ A. The proposed cutoff-independent regularization predicts a finite gap at high 
densities and indicates a smooth matching with the weak coupling QCD prediction for the gap at 
asymptotically high densities. 
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The study of the properties of high density quark mat- 
ter has attracted great interest recently - for reviews 
and extensive lists of references see Refs pj. The ear- 
lier studies 0] on this subject found energy gaps of the 
order of a few MeV's. Since gaps of this order are 
much too small to have observable consequences, not 
much attention wasjdven to the subject until recently, 
when it was shown (3, |j] within the context of instanton- 
motivated four-fermion interactions that gaps of the or- 
der of 100 MeV could be obtained. The possibility of gaps 
of this order were corroborated by subsequent study |5j 
using weak coupling renormalization group techniques for 
QCD. The result of Ref.Q for the two-flavor spin-0 su- 
perconducting gap A can be written as 



A ~ 4 exp 
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where fi is the chemical potential and g — g(fi) is the 
QCD coupling constant. This result is clearly nonper- 
turbative, but was derived assuming weak coupling, an 
assumption likely to be valid only at very high densities. 
Although inapplicable for densities typically found in the 
interiors of neutron stars, it seems to be a sound predic- 
tion for the color superconducting gap at asymptotically 
high baryon number densities. In particular, using the 
one-loop running of <?(/i) with fi, Eq. JQl predicts that 
A(/i) is an increasing function of \i. Elaborations and 
corrections Q to Eq. (JJJ do not change this behavior. 

In view of the inapplicability of weak coupling tech- 
niques at densities of phenomenological interest and the 
fact that nonperturbative lattice techniques are not yet 
sufficiently developed to deal with such problems, the 
use of phenomenological models seem to be necessary to 
make progress in the field. In this context, models with 
nonrenormalizable four-fermion interactions have been 
extensively used to study different aspects of dynami- 
cal chiral symmetry breaking (D%SB) and high density 
quark matter Jjj. However, four-fermion models at the 
one loop level predict vanishing superconducting gaps at 
high densities, a feature that is caused by the use of a 



regularizing momentum cutoff A of the divergent one- 
loop integrals [a, |y, Il0| . Since the phenomenon of super- 
conductivity involves momenta of the order of the Fermi 
momentum kp, for baryon number densities such that 
kp ~ A the cutoff regularization becomes clearly inad- 
equate and the vanishing of the gap at high densities 
might not be a physical feature of the problem. Although 
the QCD prediction of Eq. (JJJ of a nonzero gap is valid 
only at very high densities, there seems no physical mo- 
tivation for expecting vanishing gaps at moderately high 
densities. Of course, in principle one cannot exclude the 
possibility that the spin-0 gap (considered here) indeed 
vanishes at intermediate values of the chemical potential 
and then is again different from zero at much higher den- 
sities, as predicted by QCD - but such a question cannot 
be answered without a detailed nonperturbative calcula- 
tion within QCD. 

The aim of the present paper is to present an alterna- 
tive to the cutoff regularization of point-like four-fermion 
interactions. Specifically, the superconducting one-loop 
integrals are reorganized through mathematical identi- 
ties in a way finite integrals become separated from jj,- 
independent divergent integrals. The finite integrals are 
integrated without imposing any restriction to the inte- 
gration momenta and the divergent integrals are related 
to physical quantities at the D^SB scale in vacuum. As a 
result, one obtains a superconducting gap that is nonvan- 
ishing at high baryon number densities. We will keep the 
discussion very general, no specific four-fermion model 
will be used to explain our procedure; commitment to a 
specific model will be made only when presenting numer- 
ical results. 

Invariably, at the one-loop level point-like four-fermion 
interactions of massless quarks lead to a gap equation for 
the superconducting gap A of the form 



1 = XGi 



1 
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where A is a numerical coefficient that depends on the 
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Lorentz, color and flavor structure of the interaction, and 
G is the strength of the four-fermion interaction. Obvi- 
ously, the above integral is ultraviolet divergent and a 
rcgularization scheme must be specified in order to pro- 
ceed. The traditional way to regularize the integral of 
Eq. is via a three- or four-dimensional sharp cutoff, 
or via a form factor. In view of the nonrenormalizablc 
nature of the interaction, the regularization becomes part 
of the model. Using a three-dimensional sharp cutoff, the 
self-consistent solution for A = A(/i) leads to the result 
that A = for \x ~ A - the use of a form factor instead 
leads to the same qualitative behavior. 

The regularization method we advocate here, proposed 
in Ref. ^lj and used in different contexts ^3] , avoids the 
use of an explicit regulator and the calculation of any 
divergent integrals. Specifically, instead of introducing 
a cutoff in the integrals of Eq. (|2J), the integrands are 
manipulated in a way divergences are isolated in well- 
known, /i-independent one-loop divergent integrals that 
can be related to divergent integrals of the problem of 
D^SB in vacuum. Finite integrals are integrated with- 
out imposing any restriction to their integrands and the 
remaining divergent integrals fitted to physical quantities 
at the D^SB scale in vacuum. 

This is done in the following way. In order to safely ma- 
nipulate the divergent integrals, we initially assume the 
integral in Eq. @ regularized by some regulating func- 
tion /(fc 2 /A 2 ), were A is the regularization parameter. 
We rewrite the gap equation as 



1 = XGi 



d^k 



k 2 



(k + fi) — A 2 
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where J. indicates that the integral is regularized. The 
next step is the reorganization of the integrands as men- 
tioned above. This goal can be achieved by using the 
identity 

1 1 



k%-(k±n) -A 2 
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which can be obtained after using three times in succes- 
sion the identity 

1 1 



fc 2 - (fc±^) 2 - A 2 



fc 2 - fc 2 - A 2 



H 2 ± 2k[i 



(fc 2 -fc 2 -A 2 ) fc 2 -(fc±/x) 



(5) 



This identity corresponds to a subtraction of the original 
expression at the point \i = 0, akin to what is done with 



nonconvergent integrals in dispersion relations. When 
the identity of Eq. (0J is substituted into Eq. 10, one 
can rewrite the gap equation in the form 



1 = 8AG{2 [il quad (A 2 )] - V 2 [iI log (A 2 )] 
+ I fm (A 2 ,v) + / /m (A 2 ,-^)}, 



where 
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In the absence of the regularizing function, I qua d(A 2 ) and 
I; og (A 2 ) are divergent, but Ifi n (A 2 , /i) is finite. The di- 
vergent integrals are of the same form as those appearing 
in the problem of D^SB in vacuum, with the difference 
that the mass scale appearing in the superconductivity 
integrals is A while in the D^SB integrals the scale is 
the constituent quark mass M. But one can relate the 
integrals at different mass scales by making use of the 
following scaling properties [ill 
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These can be obtained by integrating the relations 

dlquad (M 2 ) _ j ( 2 \ dllogjjJ 2 ) _ ~l 

d^ 2 ~ 109 V > ' 
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between the two mass scales M 2 and A 2 . The expressions 
in Eq. \Yl\\ can be proven by direct differentiation of the 
integrals in Eqs. Q and ©. In integrating these equa- 
tions, finite integrals are integrated without regulariza- 
tion. The relations of Eqs. I|1U|) and allow us to nor- 
malize the divergent integrals I qua d (A 2 ) and Ii og (A 2 ) 
to chiral observables at the scale of the constituent quark 
mass M. This is so because I qU ad (M 2 ) and Ii og (M 2 ) in 
this class of models are related respectively to the chiral 
condensate {ipip) and to the pion decay constant f v (in 
the chiral limit) by 



lIquad(M 2 ) 



-W) 

12M : 



Ul og (M 2 ) = - 
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(13) 



Therefore, making use of the relations of Eqs. (|10|) and 
ijTTjl. one can eliminate I qua d{A 2 ) and Ii og (A 2 ) in the gap 
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equation for A in favor of (ipip) and /„., obtaining 



1 = AG 



27i(M 2 ,A 2 ) ~4fi 2 I 2 (M 2 ,A 2 



+ I fm {A 2 ,fi)+I fm {A 2 ,-fi) 
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In this way, we have eliminated the divergences with- 
out specifying any regularization function. The use of 
f(k 2 /A 2 ) is simply a matter of formality, in the sense 
that once the integrals are regularized, they can be freely 
manipulated. In a renormalizable theory, one would sim- 
ply do the same manipulations but the remaining diver- 
gences would be eliminated by counterterms added to the 
original Lagrangian (or Hamiltonian) . The basic diver- 
gent integrals could also be written in terms of an arbi- 
trary mass scale m, which can be used to relate observ- 
ables at different mass scales by using the standard renor- 
malization group techniques. The entire regularization- 
renormalization process can be done without specifying 
any regularization function and without calculating any 
divergent integral. In the present case with nonrenormal- 
izable interactions, the divergences are simply fitted to 
physical quantities at the chiral symmetry mass scale M. 
Note that instead of referring to a regularizing function 
f(k 2 /A 2 ) parameterized in terms of a momentum scale 
A, one could have simply used dimensional regulariza- 
tion, fit the divergent integrals to (^ip) and / T as above, 
and evaluate finite integrals at the physical dimension. 

So far, the discussion has been completely general, 
valid for any point-like four-fermion interaction. In or- 
der to see the physical consequences of our proposed 
method, we solve the gap equation for a specific four- 
fermion model. For illustrative purposes, we use the 
simple chirally symmetric SU (2) Nambu-Jona-Lasinio 
model with scalar-plus-pseudoscalar four-fermion inter- 
actions For such a model, A = 8Nf, where Nf = 2 
is the number of flavors. The model needs only two pa- 
rameters as input, one is the strength of the interaction 
and the other is a regularization parameter which in the 
present approach can be taken to be one of the divergent 
integrals. Using a standard value for the strength of the 
interaction, G = 3.1 GeV~ 2 , and (tpip) = (- 250 MeV) 3 
it is possible to obtain approximately M — 313 MeV 
and — 93 MeV. With these values, one obtains the 
solid line in Fig. Clearly, this results shows that the 
superconducting gap A is an increasing function of the 
chemical potential \i. Although we have not plotted the 



result for larger values of /z, the curve actually keeps in- 
creasing very slowly as a logarithmic function, and for 
/i ~ 3 GeV the gap is A ~ 300 MeV. The basic reason 
for the nonvanishing of A at high values of \i is that no 
momentum cutoff is enforced on the finite integral and 
therefore there is enough phase space at large chemical 
potentials for allowing larger gaps. To substantiate this 
we have imposed a cutoff A = 650 MeV on Ifi n and ob- 
tained the dashed curve in Fig. ^ The result is clear, 
the gap vanishes for chemical potentials of the order of 
the cutoff /i ~ A. Note that for consistency one should 
have used the same cutoff in the finite integrals that ap- 
pear in the derivation Eqs. 1)1 0|l and (| 1 1 11 - but this would 
simply worsen the situation, in the sense that A would 
be zero for fi < 700 MeV. We repeated the calculation 
using other types of point-like four-fermion interactions 
and obtained qualitatively similar results. 
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FIG. 1: Superconducting gap A as a function of the quark 
chemical potential. The solid line is the solution of Eq. (I14H 
and the dashed line is the result using a cutoff in Ia n . 

The same scheme can be used at finite tempera- 
tures ^3 j the g a P equation would be given by an ex- 
pression similar to Eq. (|14f> . with the difference that 
Ifi n (A 2 , fi) would include a term with the Fermi-Dirac 
distributions for quarks and anti-quarks. In connection 
to this, it is worth mentioning that the integration of 
finite integrals without imposing a cutoff has been ad- 
vocated in Refs. [l4| in the context of D^SB, where the 
finite-temperature integrals with the Fermi-Dirac distri- 
butions are freely integrated, while the divergent inte- 
grals are still evaluated with a cutoff A. The crucial dif- 
ference with the present approach is the use of the scaling 
relations of Eqs. (|10|) and (|11|) . by which one can always 
relate the divergent integrals to vacuum quantities with- 
out commitment to a specific regularization. Detailed fi- 
nite temperature results using the scheme proposed here 
will be presented elsewhere |l3j. 

Before concluding, one should notice that although we 
have not imposed a cutoff on finite integrals, there is still, 
of course, an implicit regularization scale dependence in 
the model through the fitting I qua d and Ii og to vacuum 
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physical quantities (t/jip) and / ff as in Eq. i|13|) . This 
dependence reflects the fact that the physical quantities 
(and the coupling G) are fitted at the implicit regular- 
ization scale A. Equivalently, had we used explicit di- 
mensional regularization, a mass scale would enter the 
problem to match the physical dimensions of the inte- 
grals and this scale would also be implicit in the fitting 
of the divergent integrals. This scale sets a limit on the 
chemical potential /i for which the fitting makes physical 
sense. Therefore, to use a four-fermion model at very 
high densities one must devise a scheme to extend the 
model beyond the vacuum scale. This actually can be 
done using the method proposed in Ref. Q by making 
the coupling run with A, i.e. by making G = G(A), and 
postulating a \x dependence of A. In Ref. || this was done 
by demanding A-independence of / w in a cutoff regular- 
ization scheme. In the present approach such a method 
can be applied using the scaling properties of the diver- 
gent integrals I qua d and Ii og . A detailed discussion on 
this will be presented in a separate publication 0|. In 
this context, it would also be interesting to implement 
the method discussed here in instanton-motivated inter- 
actions, where the cutoff has a physical origin and has a 
density dependence such that the associated form factor 
peaks around the Fermi momentum with a width of the 
order of the inverse size of the instanton |l5j . Since the 
instanton effects at large density are suppressed because 
the effective coupling becomes small, the present method 
might be useful for handling a possible smooth matching 
to perturbative QCD. 

In summary, we have shown that using a cutoff- 
independent regularization of nonrenormalizable point- 
like four-fermion interactions at the one-loop level one 
obtains nonvanishing superconducting gaps at high den- 
sities in quark matter. The result also indicates a smooth 
matching without abrupt discontinuities with the weak 
coupling QCD prediction for the gap 0, |||. Although 
this QCD prediction of a nonzero gap is valid only at 
asymptotically high quark densities, it seems neverthe- 
less reasonable to assume that the vanishing of the gap 
at high densities in four-fermion models is a nonphysical 
artifact of the cutoff regularization. Admittedly, there is 
implicit an assumption that four-fermion models can be 
used at energies higher than the chiral symmetry break- 
ing scale. This can be naturally implemented by making 
regularization parameters density dependent, as physi- 
cally motivated in Ref. In this sense, the present ap- 
proach opens new opportunities for applications of four- 
fermion models for studying properties of high density 
quark matter without facing the difficulties imposed by 
a cutoff. Obviously, the technique presented here is ap- 
plicable in a wider context where nonrenormalizable in- 
teractions are used for modelling physical systems, such 
as atomic condensates, for example. 
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